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Abstract 



We search for the decay of the meson into a pair of leptons in the sup- 
pressed channels — > e^e~, B^ — > and in the lepton number vi- 
olating channel B^ e^/x^ in a sample of 9.7 x 10^ BB pairs recorded 
by CLEO detector. No signal is found, and the following upper limits 
on the branching fractions are established: B{B^ -> e+e") < 8.3 x lO"'^, 
B{B^ fi+fi-) < 6.1 X 10"^ B{B° e^fi^) < 15 x lO""^ at 90% confidence 
level. A new lower limit on the Pati-Salam leptoquark mass Mlq > 27 TeV 
is established at 90% confidence level. 
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I. INTRODUCTION 



In the Standard Model the decays — > e^e and /x'^/i are allowed via the 

diagrams shown in Fig |l|. The Standard Model predictions for the branching fractions are 
B{B^ e+e-) = 1.9 x IQ-^^ and B{B^ = 8.0 x 10~^\ respectively [g. These 

rates put them beyond the reach of current experiments. Their observation at CLEO would 
provide clear evidence for physics beyond the Standard Model. 

In two Higgs doublet models, the corresponding branching fractions can be strongly en- 
hanced due to additional diagrams involving Higgs bosons 0. Supersymmetric particles can 
further modify the expectation, but calculations of SUSY contributions have not yet been 
made. In some models with Z-mediated flavor changing neutral currents the correspond- 
ing branching fraction can be enhanced by a factor of 400 0] compared to Standard Model 
predictions. 

The decay B^ e^/i^ is forbidden in the Standard Model by lepton number conservation. 
However it can occur in Pati-Salam leptoquark models due to the existence of particles 
that couple to both quarks and leptons. 

Currently the best 90% confidence level limits on the branching fractions are B{B^ 
e+e-) < 59 X 10"^ (from CLEO [g), B{B° eV^) < 35 x 10"^ (from CDF 0) and 
B{B^ < 6.8 X 10-^ (from CDF §). 

II. DATA SAMPLE 

We use data taken with the CLEO II/II.V detector [Q operating at the Cornell Elec- 
tron Storage Ring (CESR). The sample consists of 9.1 fb^^ taken on the T(4S') resonance, 
corresponding to approximately 9.7 x 10^ BB pairs. An additional sample of 4.5 fb~^ 
accumulated approximately 60 MeV below the T(4S') resonance is used to estimate the 
no n- resonant background. Charge conjugation is implied throughout this paper and equal 
production of charged and neutral B mesons at the T(45') resonance is assumed. The CLEO 
II charged particle tracking system consists of a 6-layer straw tube chamber, a 10-layer pre- 
cision tracking chamber and a 51 layer main drift chamber. In the CLEO II. V configuration 
the straw tube chamber was replaced by a 3-layer, double-sided silicon tracker. Beyond the 
tracking chambers, but inside the 1.5T solenoid magnet, are a time-of-fiight system and an 
electromagnetic calorimeter consisting of 7800 Csl crystals. Electron identification is per- 
formed with a likelihood function combining calorimeter {E/p, the ratio of energy deposited 
by particle in calorimeter to its momentum) and specific ionization (dE/dx) information. 
Muons are identified using proportional counters placed at various depths in the steel return 
yoke of the magnet. 

Since, at a symmetric e'^e~ collider, B mesons at the T(4S') resonance are produced 
nearly at rest, the decays B^ e~^e~, B^ f^^fJ'^ and B^ e^fi^ will contain two 
nearly back-to-back leptons with momenta p ~ 2.6 GeV. In order to select BB events, we 
require at least five charged tracks in the event with two well reconstructed and oppositely 
charged tracks identified as leptons. The lepton pair is required to have a mass within three 
standard deviations from the the mass of the B^ meson. We use the beam constrained mass, 

M{B) = ^JeI^,^ - {pi+ +pi-f. The resolution on M{B) is 2.8 MeV for the 5° e+e" 
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decay and 2.6 MeV for the 5° —>■ ^'^ ^~ and e^^"^ decays. For signal events, M{B) is 

required to be within 3 standard deviations of the 5° meson mass. 

We calculate the energy difference between the candidate meson and the beam energy 
AE = Ei+ + El- — Ebeam and require that the difference should be \AE\ < 75 MeV for 
B^ — s> fi^fi^ decay, and for 5° — >■ e"'"e~ and B^ — * e'^ff^ decays it should satisfy —100 MeV < 
AE < 75 MeV. The asymmetric AE requirement in 5° e~^e~ and B^ e^fi^ channels 
takes into account energy loss due to final state radiation from the electrons. 

The main background comes from the processes e^e~ —>■ qq{q = u, c, s, d) and e+e^ 
T^T~ . Such events usually exhibit a two jet structure and produce high momentum, ap- 
proximately back-to-back tracks which satisfy the requirements imposed on our candidate 
events. We calculate the angle 6t between the thrust axis of the B^ candidate tracks and 
the thrust axis of the remaining charged and neutral tracks in the event. The distribution of 
I cos^tI is flat for signal events and strongly peaked at 1 for background events. We require 
I cos^rl < 0.9. 



III. SIGNAL SIMULATION 



In order to take into account final state radiation (FSR) from the final state leptons, the 
matrix element for the decay J/ip e^e~(7) from Ref. |T0|, modified in order to simulate 
B^ l'^l~{pf) with /+/" = e^e~,fi'^^~ is included in the event generator. The detector 



simulation program is based on GEANT [O], the simulated events are processed as the 



data. The overall efficiencies for detecting the decays B^ —>■ e^e~ and B^ — > are 
37.3% and 46.2%, respectively, without FSR. Including FSR reduces these efficiencies to 
31.1% and 42.4%, respectively. We use the efficiencies obtained with FSR included for the 
derivation of limits on branching fractions. 

The efficiency for the decay B^ — > e^/i^ is 43.6% without FSR. There is no matrix 
element available for the decay 5° — > e^fi^lj). In this mode we use Monte Carlo without 
FSR. Since final state radiation reduces the efficiencies for B^ e~^e~ and B^ fi^fi^ 
by 20% and 9%, respectively, we estimate the effect of FSR in B^ — > e^/x^ to be 15% and 
include this in the derivation of the upper limit on the branching fraction as a systematic 
uncertainty. 



IV. RESULTS 

The background from the process e~^e~ qq results from misidentification of hadrons 
as leptons and has been determined as follows. First, the analysis is redone without impos- 
ing the lepton identification requirements to determine the total number of 2-track events 
passing all other selection criteria in the signal and sideband regions (5.260 GeV < M{B) < 
5.270 GeV and |Aii^| < 200 MeV) for both the on-resonance and off-resonance data samples. 
These are used to determine the expected number of events in the on-resonance signal region. 
(This is done by multiplying the number of events seen in the signal region in off-resonance 
data by the ratio of numbers of events seen in the sideband region in on- and off-resonance 
data.) The probabilities for these tracks to be identified as leptons is determined from the 
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relative abundances of particles (vr, proton, electron, muon) and the fake rates. The par- 
ticle abundances for the candidate tracks in qq events have been determined using Monte 
Carlo and were found to be 0.58 for vr^, 0.25 for K^, 0.15 for and 0.004 for electrons 
and muons. For vr^ and the fake rates are determined experimentally from samples of 
D*^ — * D^Ti^ — s> K^Ti^Ti^ decays. Proton fake rates are taken from Monte Carlo. The 
expected number of backgrounds events from e^e" qq are 0.01 ± 0.01, 0.14 ± 0.03 and 
0.09 ± 0.03 for the ee, /x/i and e/x channels, respectively. 

Backgrounds from e+e~ —>■ t^t~ events are estimated using a Monte Carlo sample an- 
alyzed without lepton identification, using the procedure described above to determine the 
expected rate of events with the two candidate tracks identified as leptons. The particle 
abundances for the 5° candidate tracks in simulated r+r" events were found to be 0.74 for 
n"^, 0.04 for K"^, and 0.11 for electrons and muons. The expected numbers of background 
events are 0.10 ± 0.07 , 0.08 ± 0.06 and 0.40 ± 0.20 events for the ee, /z/z and e/z channels, 
respectively. The background from BB events is estimated from Monte Carlo and is found 
to be less than 0.01 events in all modes. 

We find no events in the signal regions for the e^e~ and modes and two events in the 
signal region for the e^/x"^ mode. These events are consistent with the expected background. 
To conservatively account for the uncertainty in the efficiency determination, we reduce the 
efficiency by the estimated systematic error of 7.6% in the e+e~ and modes and 16.3% 

in the e^fi"^ mode when calculating the upper limits. The latter error is dominated by the 
uncertainty from final state radiation. A summary of the efficiencies, event yields and limits 
is given in Table |. 

TABLE I. Summary of results and branching fraction upper limits. N^^p is the number of 
expected background events, Nobs is the number of events observed, Nm is the number of events 
corresponding to a 90% CL limit, and e is the detection efficiency. The first quoted error on e is 
statistical while the second is systematic. 



mode 


Nexp 


Nobs 


NuL 


e 


Upper limit 




fiO ^ e+e" 


0.11 ±0.07 





2.30 


31.1±0.4±2.4% 


8.3 X lO"'^ 






0.22 ±0.07 





2.30 


42.4±0.5±3.2% 


6.1 X lO"'^ 






0.49 ±0.20 


2 


5.24 


43.6±0.5±7.1% 


14.9 X 10"^ 





V. LIMITS ON PATI-SALAM LEPTOQUARKS 

One of the simplest models that assumes symmetry between quarks and leptons is the 
Pati-Salam model ||^. This model predicts heavy spin-one gauge bosons, called leptoquarks 
(LQ), that carry both color and lepton quantum numbers. While most leptoquark models 
assume that the leptoquarks couple within one generation only, Pati-Salam model allows 
cross-generation couplings. Consequently, Pati-Salam leptoquarks can mediate the decay 
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FIG. 1. Feynman diagrams for the decay l^l {I = e,fi) in the Standard Model. 



e^ix^ [|T^,|l^. The relationship between Bi^B'^ e"*^//^) and leptoquark mass is flS 

1 



where 



R 



UIbo 



12 
' 23 



The first attempt to constrain the Pati-Salam leptoquark masses by measuring the 
branching fraction for B^ — > e^fi^ was made by the CDF collaboration 0. The result 
obtained was Mlq > 21.7 TeV. We follow the CDF procedure by using Fb = 175 ± 30 MeV 
for 5° decay constant H, tubo = 5279.1 ±0.7(stat.) ±0.3(syst.) MeV [|l|], rrit = 176.0±6.5 
GeV [|T5| and obtain a new limit on the Pati-Salam leptoquark mass of Mlq > 27 TeV at 
90% confidence level. 

In conclusion we find no evidence for the decays B^ e^e~ , fi^fi^ and e^/x^ and obtain 
upper limits on the corresponding branching fractions of B{B^ e^e~) < 8.3 x 10^'', 
B{B^ < 6.1 X 10"^ B{B^ e±/i^) < 15 x 10"^ at 90% confidence level. The 

limit on B{B'^ —>■ e'^fi^) corresponds to a lower limit of 27 TeV on the Pati-Salam leptoquark 
mass at 90% CL. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. I.P.J. Shipsey thanks the NYI program of the NSF, 
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FIG. 2. Beam constrained mass of dilepton pair vs energy difference AE = Ei+ + — Eh^am 
for on resonance data. All selection criteria except the requirements on lepton pair mass and AE 
are included. The box in the center of each plot shows the signal region. 
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